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Abstract — It is well-known that asymptotic stability (AS) of 
homogeneous polynomial vector fields of degree one (i.e., linear 
systems) can be decided in polynomial time e.g. by searching for 
a quadratic Lyapunov function. Since homogeneous vector fields 
of even degree can never be AS, the next interesting degree to 
consider is equal to three. In this paper, we prove that deciding 
AS of homogeneous cubic vector fields is strongly NP-hard and 
pose the question of determining whether it is even decidable. As 
a byproduct of the reduction that establishes our NP-hardness 
result, we obtain a Lyapunov-inspired technique for proving 
positivity of forms. We also show that for asymptotically stable 
homogeneous cubic vector fields in as few as two variables, the 
minimum degree of a polynomial Lyapunov function can be 
arbitrarily large. Finally, we show that there is no monotonicity 
in the degree of polynomial Lyapunov functions that prove AS; 
i.e., a homogeneous cubic vector field with no homogeneous 
polynomial Lyapunov function of some degree d can very well 
have a homogeneous polynomial Lyapunov function of degree 
less than d. 

I. Introduction 

A. Background 

We are concerned in this paper with a continuous time 
dynamical system 

i = f{x), (1) 

where / : M" M" is a polynomial and has an equilibrium 
at the origin, i.e., /(O) = 0. Polynomial differential equa- 
tions appear ubiquitously in engineering and sciences either 
as true models of physical systems, or as approximations 
to other families of nonlinear dynamics. The problem of 
deciding stability of equilibrium points of such systems is of 
fundamental importance in control theory. The goal of this 
paper is to demonstrate some of the difficulties associated 
with answering stability questions about polynomial vector 
fields in terms of both computational complexity and non- 
existence of "simple" Lyapunov functions, even if one limits 
attention to very restricted settings. 

The notion of stability of interest in this paper is (local 
or global) asymptotic stability. The origin of ([TJ is said to 
be stable in the sense of Lyapunov if for every e > 0, there 
exists a (5 = (5(e) > such that 

||a;(0)|| <5 ^ \\x{t)\\ < e, > 0. 

We say that the origin is asymptotically stable (AS) if it is 
stable in the sense of Lyapunov and 5 can be chosen such 
that 



< 5 



lim x{t) = 0. 
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The origin is globally asymptotically stable (GAS) if it 
is stable in the sense of Lyapunov and Va:(0) G R", 
limi^oo x{t) = 0. 

The degree of the vector field in ([T]i is defined to be the 
largest degree of the components of /. Our focus in this 
paper is on homogeneous polynomial vector fields. A scalar 
valued function p : R" — > R is said to be homogeneous 
(of degree d) if it satisfies p{\x) = X'^p{x) for all x G R" 
and all A G R. A homogeneous polynomial is also called a 
form. All monomials of a form share the same degree. We 
say that the vector field / in ([T]! is homogeneous if all com- 
ponents of / are forms of the same degree. Homogeneous 
systems are extensively studied in the literature on nonlinear 
control; see e.g. [1], [2], [3], [4], [5], [6], [7]. Since our 
results are negative in nature, their validity for homogeneous 
polynomial systems obviously also implies their validity for 
all polynomial systems. 

A basic fact about homogeneous vector fields is that for 
these systems the notions of local and global asymptotic 
stability are equivalent. Indeed, a homogeneous vector field 
of degree d satisfies f{Xx) ~ X'^f{x) for any scalar A, and 
therefore the value of / on the unit sphere determines its 
value everywhere. It is also well-known that an asymptot- 
ically stable homogeneous system admits a homogeneous 
Lyapunov function [8, Sec. 57], [6]. 

B. An open question of Arnold 

It is natural to ask whether stability of equilibrium points 
of polynomial vector fields can be decided in finite time. In 
fact, this is a well-known question of Arnold that appears 
in [91: 

'Is the stability problem for stationary points algorith- 
mically decidable? The well-known Lyapunov theorerrQ 
solves the problem in the absence of eigenvalues with 
zero real parts. In more complicated cases, where the 
stability depends on higher order terms in the Taylor 
series, there exists no algebraic criterion. 
Let a vector field be given by polynomials of a fixed 
degree, with rational coefficients. Does an algorithm 
exist, allowing to decide, whether the stationary point 
is stable?" 

To our knowledge, there has been no formal resolution to 
this question, neither for the case of stability in the sense 
of Lyapunov, nor for the case of asymptotic stability (in 
its local or global version). In [10], da Costa and Doria 

'The theorem that Arnold is referring to here is the indirect method 
of Lyapunov related to linearization. This is not to be confused with 
Lyapunov's direct method (or the second method), which is what we are 
concerned with in sections that follow. 



show that if the right hand side of the differential equation 
contains elementary functions (sines, cosines, exponentials, 
absolute value function, etc.), then there is no algorithm for 
deciding whether the origin is stable or unstable. They also 
present a dynamical system in [11] where one cannot decide 
whether a Hopf bifurcation will occur or whether there will 
be parameter values such that a stable fixed point becomes 
unstable. A relatively larger number of undecidability re- 
sults are available for questions related to other properties 
of polynomial vector fields, such as reachability [12] or 
boundedness of domain of definition [13], or for questions 
about stability of hybrid systems [14], [15], [16], [17]. We 
refer the interested reader to the survey papers in [18], [12], 
[19], [20], [21]. 

We are also interested to know whether the problem of 
deciding asymptotic stability of homogeneous polynomial 
vector fields is undecidable for some fixed degree, say, 
equal to 3. The answer to such decidability questions, or 
at least the level of difficulty associated with proving such 
results, can depend in a subtle way on the exact criteria in 
question. For example, it has been known for a while that 
the question of determining boundedness of trajectories for 
arbitrarily switched linear systems is undecidable [15] even 
when one restricts attention to switched systems defined by 
nonnegative matrices. On the other hand, the complexity of 
testing asymptotic stability for the same class of systems 
remains open and in fact is conjectured to be decidable [22]. 

C. Existence of polynomial Lyapunov functions 

For stability analysis of polynomial vector fields, it is 
most common (and quite natural) to search for Lyapunov 
functions that are polynomials themselves. This approach has 
become further prevalent over the past decade due to the fact 
that techniques from sum of squares optimization [23] have 
provided for algorithms that given a polynomial system can 
efficiently search for a polynomial Lyapunov function [23], 
[24]. The question is therefore naturally motivated to de- 
termine whether stable polynomial systems always admit 
polynomial Lyapunov functions, and whether one can give 
upper bounds on the degree of such Lyapunov functions in 
cases when they do exist. A study of questions of this type 
for different notions of stability has recently been carried 
out in [25], [26], [27], [28], [29, Chap. 4]. In this paper, we 
continue this line of research by studying the case where the 
vector field is homogeneous. 

Throughout this paper, by a (polynomial) Lyapunov func- 
tion for ([T]i, we mean a positive definite polynomial function 
V whose derivatives V along trajectories of ([U is negative 
definite; i.e., a function V satisfying 

V{x) > V.T7^0 (2) 
V{x) = {VV{x)J{x)) < Vs^O. (3) 

Here, W{x) denotes the gradient vector of V, and (.,.) 
is the standard inner product in R". If such a y is also 
radially unbounded, then the inequalities in (|2]i and (O imply 
that the origin of ([T]i is GAS. When the dynamics / is 
homogeneous, we can restrict our search to homogeneous 



polynomials. Such a Lyapunov function is automatically 
radially unbounded and proves (local or equivalently global) 
asymptotic stability of the homogeneous vector field. 

Naturally, questions regarding complexity of deciding 
asymptotic stability and questions about existence of Lya- 
punov functions are related. For instance, if one proves that 
for a class of polynomial vector fields, asymptotic stability 
implies existence of a polynomial Lyapunov function to- 
gether with a computable upper bound on its degree, then 
the question of asymptotic stability for that class becomes 
decidable. This is due to the fact that given a polynomial 
system and an integer d, the question of deciding whether 
the system admits a polynomial Lyapunov function of degree 
d can be answered in finite time using quantifier elimina- 
tion [30], [31]. 

For the case of linear systems (i.e., homogeneous systems 
of degree 1), the situation is particularly nice. If such a 
system is asymptotically stable, then there always exists 
a quadratic Lyapunov function. Asymptotic stability of a 
linear system x — Ax is equivalent to the easily checkable 
algebraic criterion that the eigenvalues of A be in the open 
left half complex plane. Deciding this property of the matrix 
A can formally be done in polynomial time, e.g. by solving 
a Lyapunov equation [21]. 

Moving up in the degree, it is not difficult to show that 
if a homogeneous polynomial vector field has even degree, 
then it can never be asymptotically stable; see e.g. [8, p. 
283]. So the next interesting case occurs for homogeneous 
vector fields of degree 3. We will prove three results in this 
paper which demonstrate that already for cubic homogeneous 
systems, the situation is significantly more complex than it 
is for linear systems. We outline our contributions next. 

D. Contributions and organization of this paper 

In Section nil we prove that determining asymptotic stabil- 
ity for homogeneous cubic vector fields is strongly NP-hard 
(Theorem 12. II ). Although this of course does not resolve the 
question of Arnold, the result gives a lower bound on the 
complexity of this problem. It is an interesting open question 
to investigate whether in this specific setting, the problem is 
also undecidable. 

The implication of the NP-hardness of this problem is that 
unless P=NP, it is impossible to design an algorithm that can 
take as input the (rational) coefficients of a homogeneous cu- 
bic vector field, have running time bounded by a polynomial 
in the number of bits needed to represent these coefficients, 
and always output the correct yes/no answer on asymptotic 
stability. Moreover, the fact that our NP-hardness result is in 
the strong sense (as opposed to weakly NP-hard problems 
such as KNAPSACK, SUBSET SUM, etc.) impHes that the 
problem remains NP-hard even if the size (bit length) of the 
coefficients is 0(log7i), where n is the dimension. For a 
strongly NP-hard problem, even a pseudo-polynomial time 
algorithm cannot exist unless P=NP. See [32] for precise 
definitions and more details. 

In Section [III we also present a Lyapunov-inspired tech- 
nique for proving positivity of forms that comes directly 



out of the reduction in the proof of our NP-hardness result 
(Corollary 12. 11 1. We show the potential advantages of this 
technique over standard sum of squares techniques on an 
example (Example 12. 11 1. 

In Section Uni we prove that unlike AS linear systems that 
always admit quadratic Lyapunov functions, AS cubic homo- 
geneous systems may need polynomial Lyapunov functions 
of arbitrarily large degree, even when the dimension is fixed 
to 2 (Theorem l3.1l ). Finally, in Section llVl we show that there 
is no monotonicity in the degree of homogeneous polynomial 
Lyapunov functions for homogeneous cubic vector fields. 
We give an example of such a vector field which admits 
a homogeneous polynomial Lyapunov function of degree 4 
but not one of degree 6 (Theorem 14. It . 

II. NP-HARDNESS OF DECIDING ASYMPTOTIC STABILITY 
OF HOMOGENEOUS CUBIC VECTOR FIELDS 

The main result of this section is the following theorem. 

Theorem 2.1: Deciding asymptotic stability of homoge- 
neous cubic polynomial vector fields is strongly NP-hard. 

The key idea behind the proof of this theorem is the 
following: We will relate the solution of a combinatorial 
problem not to the behavior of the trajectories of a cubic 
vector field that are hard to get a handle on, but instead 
to properties of a Lyapunov function that proves asymptotic 
stability of this vector field. As we will see shortly, insights 
from Lyapunov theory make the proof of this theorem quite 
simple. The reduction is broken into two steps: 

ONE-IN-THREE 3SAT 

; 

positivity of quartic forms 

; 

asymptotic stability of cubic vector fields 

A. Reduction from ONE-IN-THREE 3 SAT to positivity of 
quartic forms 

A form q is said to be nonnegative or positive semidefinite 
if q{x) > for all x in M". We say that a form q is positive 
definite if q{x) > for all a; 7^ in R". (Note that forms 
necessarily vanish at the origin.) It is well-known that decid- 
ing nonnegativity of quartic forms is NP-hard; see e.g. [33] 
and [34]. For reasons that will become clear shortly, we are 
interested instead in showing hardness of deciding positive 
definiteness of quartic forms. This is in some sense even 
easier to accomplish. A very straightforward reduction from 
3SAT proves NP-hardness of deciding positive definiteness 
of polynomials of degree 6. By using ONE-IN-THREE 3SAT 
instead, we will reduce the degree of the polynomial from 6 
to 4. 

Proposition 1: It is strongljH NP-hard to decide whether 
a homogeneous polynomial of degree 4 is positive definite. 

Proof: We give a reduction from ONE-IN-THREE 
3SAT which is known to be NP-complete [32, p. 259]. Recall 
that in ONE-IN-THREE 3SAT, we are given a 3SAT instance 

-The NP-hardness results of this section will all be in the strong sense. 
From here on, we drop the prefix "strong" for brevity. 



(i.e., a collection of clauses, where each clause consists of 
exactly three literals, and each literal is either a variable or its 
negation) and we are asked to decide whether there exists a 
{0, 1} assignment to the variables that makes the expression 
true with the additional property that each clause has exactly 
one true literal. 

To avoid introducing unnecessary notation, we present the 
reduction on a specific instance. The pattern will make it 
obvious that the general construction is no different. Given 
an instance of ONE-IN-THREE 3SAT, such as the following 

(xi V X2 V X4 ) A (x2 V X3 V 2:5 ) A (xi V .T3 V .T5 ) A (a;i V 2:3 V X4 ) , 

(4) 

we define the quartic polynomial p as follows: 

+ {X1 + (1 - X2) +X4- 1)2 + ((1 - X2) 

+ - X3) + xr, - 1)^ 

+ ((l-a;i) + a;3 + (1-2:5) -1)2 

+ (2:1 + 2:3 + 2:4 - 1)2. 

(5) 

Having done so, our claim is that ^(2;) > for all 2; € 
(or generally for all x e M") if and only if the ONE-IN- 
THREE 3SAT instance is not satisfiable. Note that p is a 
sum of squares and therefore nonnegative. The only possible 
locations for zeros of p are by construction among the points 
in {0, 1}^. If there is a satisfying Boolean assignment 2; 
to (HI with exactly one true literal per clause, then p will 
vanish at point x. Conversely, if there are no such satisfying 
assignments, then for any point in {0, 1}^, at least one of 
the terms in (|5]i will be positive and hence p will have no 
zeros. 

It remains to make p homogeneous. This can be done via 
introducing a new scalar variable y. If we let 

Phix,y) = y^pi^), (6) 

then we claim that ph (which is a quartic form) is positive 
definite if and only if p constructed as in (|5]l has no zeros0 
Indeed, if p has a zero at a point x, then that zero is inherited 
by ph at the point (2;, 1). If p has no zeros, then ^ shows 
that Ph can only possibly have zeros at points with y = 0. 
However, from the structure of p in (|5]) we see that 

Ph{x,0) = 2;f H 1-2;|, 

which cannot be zero (except at the origin). This concludes 
the proof. ■ 

B. Reduction from positivity of quartic forms to asymptotic 
stability of cubic vector fields 

We now present the second step of the reduction and finish 
the proof of Theorem 12.11 

^^In general, the homogenization operation in (6) does not preserve 
positivity. For example, as shown in [35], the polynomial + (1 — 2:12:2)^ 
has no zeros, but its homogenization x'^y^ + (y^ — xiX2)^ has zeros at 
the points (1,0,0)^ and (0,1,0)^. Nevertheless, positivity is preserved 
under homogenization for the special class of polynomials constructed in 
this reduction, essentially because polynomials of type ^5) have no zeros at 
infinity. 



Proof: [Proof of Theorem 12.11 We give a reduction 
from the problem of deciding positive definiteness of quartic 
forms, whose NP-hardness was established in Proposition [T] 
Given a quartic form V V{x), we define the polynomial 
vector field 

X ^ -Vy(x). (7) 

Note that the vector field is homogeneous of degree 3. We 
claim that the above vector field is (locally or equivalently 
globally) asymptotically stable if and only if V is positive 
definite. First, we observe that by construction 

Vix) = {\7V{x),x) = -\\VV{x)f < 0. (8) 

Suppose V is positive definite. By Euler's identity for ho- 
mogeneous functions^ we have V{x) = jx'^yV{x). There- 
fore, positive definiteness of V implies that WV{x) cannot 
vanish anywhere except at the origin. Hence, V{x) < for 
all x ^ 0. In view of Lyapunov's theorem (see e.g. [36, p. 
124]), and the already mentioned fact that a positive definite 
homogeneous function is radially unbounded, it follows that 
the system in (|7]i is globally asymptotically stable. 

For the converse direction, suppose (|7]i is GAS. Our 
first claim is that global asymptotic stability together with 
V{x) < implies that V must be positive semidefinite. 
This follows from the following simple argument, which we 
have also previously presented in [37] for a different purpose. 
Suppose for the sake of contradiction that for some x E R" 
and some e > 0, we had V{x) = — e < 0. Consider a 
trajectory x{t;x) of system ^ that starts at initial condition 
X, and let us evaluate the function V on this trajectory. Since 
V{x) = — e and V{x) < 0, we have V{x{t;x)) < — e for 
all t > 0. However, this contradicts the fact that by global 
asymptotic stability, the trajectory must go to the origin, 
where V, being a form, vanishes. 

To prove that V is positive definite, suppose by con- 
tradiction that for some nonzero point x* £ R" we had 
V{x*) ~ 0. Since we just proved that V has to be positive 
semidefinite, the point x* must be a global minimum of V. 
Therefore, as a necessary condition of optimality, we should 
have W{x*) ~ 0. But this contradicts the system in Q 
being GAS, since the trajectory starting at x* stays there 
forever and can never go to the origin. ■ 

Perhaps of independent interest, the reduction we just 
gave suggests a method for proving positive definiteness 
of forms. Given a form V, we can construct a dynamical 
system as in and then any method that we may have 
for proving stability of vector fields (e.g. the use of various 
kinds of Lyapunov functions) can serve as an algorithm for 
proving positivity of V. In particular, if we use a polynomial 
Lyapunov function W to prove stability of the system in 
we get the following corollary. 

Corollary 2.1: Let V and W be two forms of possibly 
different degree. If W is positive definite, and {VW, W) is 
positive definite, then V is positive definite. 

''Euler's identity is easily derived by differentiating both sides of the 
equation V{\x) = A'*y(x) with respect to A and setting A = 1. 



A polynomial p is said to be a sutn of squares (sos) if it 
can be written as p = qf for some polynomials qi. An 

SOS polynomial is clearly nonnegative. Moreover, unlike the 
property of nonnegativity that is NP-hard to check, existence 
of an SOS decomposition can be cast as a semidefinite 
program [38], which can be solved efficiently. However, not 
every nonnegative polynomial is a sum of squares. 

An interesting fact about Corollary [2T| is that its algebraic 
version with sum of squares replaced for positivity is not true. 
In other words, we can have W sos (and positive definite), 
(VVF, W) SOS (and positive definite), but V not sos. This 
gives us a way of proving positivity of some polynomials that 
are not sos, using only sos certificates. Given a form V, since 
the expression {VW, W) is linear in the coefficients of W, 
we can use semidefinite programming to search for a form 
W that satisfies W sos and {VW, W) sos, and this would 
prove positivity of V. The following example demonstrates 
the potential usefulness of this approach. 

Example 2.1: Consider the following form of degree 6: 

V{x) — x\x\ + xfx2 — 'ix\x\x\ + x\ + TT77;(^1 +^2 + ^3)^- 

(9) 

One can check that this polynomial is not a sum of squares. 
(In fact, this is the celebrated Motzkin form [39] slightly per- 
turbed.) On the other hand, we can use the software package 
YALMIP [40] together with the SDP solver SeDuMi [41] to 
search for a form W satisfying 

W SOS 

(viy,vy) SOS. ^ ' 

If we parameterize as a quadratic form, no feasible 
solution will be returned form the solver. However, when 
we increase the degree of W from 2 to 4, the solver returns 
the following polynomial 

W{x) = 9x* + 9xf - ^x\x\ + %x\x\ + %x\x\ + 3a;| 

— x\x2 — X\x\ — x\xj, — Zx\x2X-i — 'iX\x\x-i 

— x\x'i — 4X1X2X3 — X\x\ — X2X3 

that satisfies both sos constraints in ( fTOl i. One can easily 
infer from the sos decompositions (e.g. by checking positive 
definiteness of the associated "Gram matrices") that the 
forms W and (VW, Vy) are positive definite. Hence, by 
Corollary 12.11 we have a proof that F in (|9j is positive 
definite. A 

Interestingly, approaches of this type that use gradient 
information for proving positivity of polynomials with sum 
of squares techniques have been studied by Nie, Demmel, 
and Sturmfels in [42], though the derivation there is not 
inspired by Lyapunov theory. 

111. Non-existence of a uniform bound on the 
degree of polynomial lyapunov functions in 
fixed dimension and degree 

For polynomial vector fields in general, existence of a 
polynomial Lyapunov function is not necessary for global 
asymptotic stability. In joint work with M. Krstic and PA. 
Parrilo [27], we recently gave a remarkably simple example 



of a (non-homogeneous) quadratic polynomial vector field in 
two variables that is GAS but does not admit a polynomial 
Lyapunov function (of any degree). An independent earlier 
example that appears in a book by Bacciotti and Rosier [43, 
Prop. 5.2] was brought to our attention after our work was 
submitted. We refer the reader to [27] for a discussion on 
the differences between the two examples, the main one 
being that the example in [43] does not admit a polynomial 
Lyapunov function even locally but unlike the example 
in [27] relies on using irrational coefficients. 

The situation for homogeneous polynomial vector fields, 
however, seems to be different. We conjecture that for such 
systems, existence of a homogeneous polynomial Lyapunov 
function is necessary and sufficient for (global) asymptotic 
stability. The reason for this conjecture is that we expect that 
one should be able to approximate a continuously differen- 
tiable Lyapunov function with a polynomial one on the unit 
sphere, which by homogeneity should be enough to imply 
the Lyapunov inequalities everywhere. A formal treatment 
of this idea is left for future work. Here, we build on the 
result in [43, Prop. 5.2] to prove that the minimum degree 
of a polynomial Lyapunov function for an AS homogeneous 
vector field can be arbitrarily large even when the degree and 
dimension are fixed respectively to 3 and 2. 

Proposition 2 ( [43, Prop. 5.2-a]): Consider the vector 
field 



y = 4Xx{x'^ + y'^) + 2x(2x'^ + y'^) 



(11) 



parameterized by the scalar A > 0. For all values of A the 
origin is a center for ( fTTT ). but for any irrational value of A 
there exists no polynomial function V satisfying V{x,y) = 



Theorem 3.1: Let A be a positive irrational real number 
and consider the following homogeneous cubic vector field 
parameterized by the scalar 6: 



x\ _ /cos(e) -sin(0) 
ij) ~ \sin(e) cos(0) 



-2A3/(x2+3/2)_2y(2x2+y2) 
4Ax(x2 + y2) _|_ 2a;(2x2 + y2j 



(12) 

Then for any even degree d of a candidate polynomial 
Lyapunov function, there exits a > small enough such 
that the vector field in (fT2] i is asymptotically stable but does 
not admit a polynomial Lyapunov function of degree < d. 

Proof: Consider the (non-polynomial) positive definite 
Lyapunov function 

V{x,y)^{2x^+y^t{x^ + y^) 

whose derivative along the trajectories of (fT2]) is equal to 

V{x,y) = -sin(0)(2a-2 +y2)A-i(i2 ^^2)^ 

Since V is negative definite for Q < 6 < -k, \i follows that 
for 6 in this range, the origin of (fT2l i is asymptotically stable. 

To establish the claim in the theorem, suppose for the 
sake of contradiction that there exists an upper bound d such 
that for all < 6' < TT the system admits a (homogeneous) 
polynomial Lyapunov function of degree d{6) with d{0) < d. 
Let d be the least common multiplier of the degrees d{d) 



for < 9 < TT. (Note that d{d) can at most range over 
all even positive integers less than or equal to d.) Since 
positive powers of Lyapunov functions are valid Lyapunov 
functions, it follows that for every < 9 < n, the system 
admits a homogeneous polynomial Lyapunov function Wg 
of degree d. By rescaling, we can assume without loss of 
generality that all Lyapunov functions Wg have unit area on 
the unit sphere. Let us now consider the sequence {Wg} as 
6* — 0. We think of this sequence as residing in a compact 

subset of MV d ) associated with the set j of (coefficients 
of) all nonnegative bivariate homogeneous polynomials of 
degree d with unit area on the unit sphere. Since every 
bounded sequence has a converging subsequence, it follows 
that there must exist a subsequence of {Wg} that converges 
(in the coefficient sense) to some polynomial Wq belonging 
to J. Since convergence of this subsequence also implies 
convergence of the associated gradient vectors, we get that 



Woix,y) 



dWo . ^dWo. 

-X H — - — y < 0. 



dx 



dy 



On the other hand, when 9 = 0, the vector field in (fT2b 
is the same as the one in (fTTT i and hence the trajectories 
starting from any nonzero initial condition go on periodic 
orbits. This however implies that W = everywhere and in 
view of Proposition |2] we have a contradiction. ■ 
Remark 3.1: Unlike the result in [43, Prop. 5.2], it is easy 
to establish the result of Theorem 13.11 without having to use 
irrational coefficients in the vector field. One approach is 
to take an irrational number, e.g. tt, and then think of a 
sequence of vector fields given by (fT2] i that is parameterized 
by both 9 and A. We let the fc-th vector field in the 
sequence have f?^ = \ and A^ equal to a rational number 
representing tt up to fc decimal digits. Since in the limit 
as fc ^ c« we have 0^ — > and Afe tt, it should be 
clear from the proof of Theorem 13.11 that for any integer d, 
there exists an AS bivariate homogeneous cubic vector field 
with rational coefficients that does not have a polynomial 
Lyapunov function of degree less than d. 

IV. Lack of monotonicity in the degree of 

POLYNOMIAL LYAPUNOV FUNCTIONS 

If a dynamical system admits a quadratic Lyapunov func- 
tion V, then it clearly also admits a polynomial Lyapunov 
function of any higher even degree (e.g. simply given by 
V'^ for k = 2,3,.. .). However, our next theorem shows 
that for homogeneous systems that do not admit a quadratic 
Lyapunov function, such a monotonicity property in the 
degree of polynomial Lyapunov functions may not hold. 

Theorem 4.1: Consider the following homogeneous cubic 
vector field parameterized by the scalar 9: 



-sm{9) cos{9)\ fx 



cos(9) — sm{9) J \y 



(13) 



There exists a range of values for the parameter 9 > 
for which the vector field is asymptotically stable, has no 
homogeneous polynomial Lyapunov function of degree 6, 



but admits a homogeneous polynomial Lyapunov function 
of degree 4. 

Proof: Consider the positive definite Lyapunov function 

V{x,y)=x^ + y\ (14) 
The derivative of this Lyapunov function is given by 

V{x,y) ^ -4siii(6')(a;6 

which is negative definite for < 6 < tt. Therefore, when 
6 belongs to this range, the origin of (fT2] i is asymptotically 
stable and the system admits the degree 4 Lyapunov function 
given in (fT4l i. On the other hand, we claim that for 9 small 
enough, the system cannot admit a degree 6 (homogeneous) 
polynomial Lyapunov function. To argue by contradiction, 
we suppose that for arbitrarily small and positive values of 9 
the system admits sextic Lyapunov functions Wg. Since the 
vector field satisfies the symmetry 

fx{y, -x)\ ^ / 1\ fx\ 
\y{y,~x)) \-l O) {yj ' 

we can assume that the Lyapunov functions We satisfy the 
symmetry Wg{y, —x) ~ W0{x, y). El This means that Wg can 
be parameterized with no odd monomials, i.e., in the form 

Weix, y) = cix^ + C2X^y^ + c^x'^y'^ + c^y^ , 

where it is understood that the coefficients ci , . . . , C4 are 
a function of 9. Since by our assumption Wg is negative 
definite for 9 arbitrarily small, an argument identical to the 
one used in the proof of Theorem 13.11 implies that as — !> 0, 
Wg converges to a nonzero sextic homogeneous polynomial 
Wb whose derivative Wq along the trajectories of (fTsT l (with 
= 0) is non-positive. However, note that when 9 ~ Q, the 
trajectories of (fTsT i go on periodic orbits tracing the level sets 
of the function x^ + y^. This implies that Wo = ^^V^ + 
^^-{—x"^) — 0. If we write out this equation, we obtain 

Wq = {&C\—4:C2)x'^y'^ + 2c2xy''~2c:ix'y + {Acz — (jC4)x'^y'^ — 0, 

which implies that ci = C2 = C3 = C4 = 0, hence a 
contradiction. ■ 
Remark 4.1: We have numerically computed the range 
< 6* < 0.0267, for which the conclusion of Theorem 14. 1 1 
holds. This bound has been computed via sum of squares 
relaxation and semidefinite programming (SDP) by using 
the SDP solver SeDuMi [41]. What allows the search for a 
Lyapunov function for the vector field in ( fT3] l to be exactly 
cast as a semidefinite program is the fact that all nonnegative 
bivariate forms are sums of squares. 
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^To see this, note that any Lyapunov function Vg for this system can be 
made into one satisfying this symmetry by letting We{x, y) = Ve{x,y) + 
Vgiy, -x) + Vg{-x, -y) + Vg{-y,x). 
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